A semi-analytical solution for simulating contaminant transport subject to chain-decay reactions.

We present a set of new, semi-analytical solutions to simulate three-dimensional contaminant transport subject to first-order chain-decay reactions. The aquifer is assumed to be areally semi-infinite, but finite in thickness. The analytical solution can treat the transformation of contaminants into daughter products, leading to decay chains consisting of multiple contaminant species and various reaction pathways. The solution in its current form is capable of accounting for up to seven species and four decay levels. The complex pathways are represented by means of first-order decay and production terms, while branching ratios account for decay stoichiometry. Besides advection, dispersion, bio-chemical or radioactive decay and daughter product formation, the model also accounts for sorption of contaminants on the aquifer solid phase with each species having a different retardation factor. First-type contaminant boundary conditions are utilized at the source (x=0 m) and can be either constant-in-time for each species, or the concentration can be allowed to undergo first-order decay. The solutions are obtained by exponential Fourier, Fourier cosine and Laplace transforms. Limiting forms of the solutions can be obtained in closed form, but we evaluate the general solutions by numerically inverting the analytical solutions in exponential Fourier and Laplace transform spaces. Various cases are generated and the solutions are verified against the HydroGeoSphere numerical model.

Incorporating layer- and local-scale heterogeneities in numerical simulation of unsaturated flow and tracer transport.

This study characterizes layer- and local-scale heterogeneities in hydraulic parameters (i.e., matrix permeability and porosity) and investigates the relative effect of layer- and local-scale heterogeneities on the uncertainty assessment of unsaturated flow and tracer transport in the unsaturated zone of Yucca Mountain, USA. The layer-scale heterogeneity is specific to hydrogeologic layers with layerwise properties, while the local-scale heterogeneity refers to the spatial variation of hydraulic properties within a layer. A Monte Carlo method is used to estimate mean, variance, and 5th, and 95th percentiles for the quantities of interest (e.g., matrix saturation and normalized cumulative mass arrival). Model simulations of unsaturated flow are evaluated by comparing the simulated and observed matrix saturations. Local-scale heterogeneity is examined by comparing the results of this study with those of the previous study that only considers layer-scale heterogeneity. We find that local-scale heterogeneity significantly increases predictive uncertainty in the percolation fluxes and tracer plumes, whereas the mean predictions are only slightly affected by the local-scale heterogeneity. The mean travel time of the conservative and reactive tracers to the water table in the early stage increases significantly due to the local-scale heterogeneity, while the influence of local-scale heterogeneity on travel time gradually decreases over time. Layer-scale heterogeneity is more important than local-scale heterogeneity for simulating overall tracer travel time, suggesting that it would be more cost-effective to reduce the layer-scale parameter uncertainty in order to reduce predictive uncertainty in tracer transport.

Live/real time three-dimensional transthoracic echocardiographic assessment of left ventricular volumes, ejection fraction, and mass compared with magnetic resonance imaging.

Due to reliance upon geometric assumptions and foreshortening issues, the traditionally utilized transthoracic two-dimensional echocardiography (2DTTE) has shown limitations in assessing left ventricular (LV) volume, mass, and function. Cardiac magnetic resonance imaging (MRI) has shown potential in accurately defining these LV characteristics. Recently, the emergence of live/real time three-dimensional (3D) TTE has demonstrated incremental value over 2DTTE and comparable value with MRI in assessing LV parameters. Here we report 58 consecutive patients with diverse cardiac disorders and clinical characteristics, referred for clinical MRI studies, who were evaluated by cardiac MRI and 3DTTE. Our results show good correlation between the two modalities.

A mountain-scale thermal-hydrologic model for simulating fluid flow and heat transfer in unsaturated fractured rock.

A multidimensional, mountain-scale, thermal-hydrologic (TH) numerical model is presented for investigating unsaturated flow behavior in response to decay heat from the proposed radioactive waste repository in the Yucca Mountain unsaturated zone (UZ), The model, consisting of both two-dimensional (2-D) and three-dimensional (3-D) representations of the UZ repository system, is based on the current repository design, drift layout, thermal loading scenario, and estimated current and future climate conditions. This mountain-scale TH model evaluates the coupled TH processes related to mountain-scale UZ flow. It also simulates the impact of radioactive waste heat release on the natural hydrogeological system, including heat-driven processes occurring near and far away from the emplacement tunnels or drifts. The model simulates predict thermally perturbed liquid saturation, gas- and liquid-phase fluxes, and water and rock temperature elevations, as well as the changes in water flux driven by evaporation/condensation processes and drainage between drifts. These simulations provide insights into mountain-scale thermally perturbed flow fields under thermal loading conditions.

A triple-continuum approach for modeling flow and transport processes in fractured rock.

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain.

Development of discrete flow paths in unsaturated fractures at Yucca Mountain.

We have carried out numerical modeling studies to investigate the development of discrete-fracture flow paths and flow-focusing phenomena in the unsaturated rock of the potential repository horizon at Yucca Mountain, Nevada. These studies are based on two- and three-dimensional (2-D and 3-D) numerical models using site-specific parameters. The 2-D and 3-D models use high-resolution spatial discretization to explicitly include effects of discrete fractures with stochastically developed fracture permeabilities and a continuum approach. The permeability field is generated based on air permeability measurements at various scales. For most of the cases considered, uniform infiltration with different average rates (1-500 mm/year) is prescribed at the top of the model, while variability in outflow at the bottom of the model is used to evaluate the degree of flow focusing. In addition, scenarios involving nonuniform infiltration at the top boundary, different permeability correlation lengths and different flow-allocation schemes were analyzed. The modeling results obtained from all of the cases showed a remarkably similar flow-focusing pattern at the repository horizon. Furthermore, tracer transport simulation results also revealed additional features of focused flow and transport through the fracture network.

Analysis of flow behavior in fractured lithophysal reservoirs.

This study develops a mathematical model for the analysis of pressure behavior in fractured lithophysal reservoirs. The lithophysal rock is described as a tri-continuum medium, consisting of fractures, rock matrices, and cavities. In the conceptual model, fractures have homogeneous properties throughout and interact with rock matrices and cavities that have different permeabilities and porosities. Global flow occurs through the fracture network only, while rock matrices and cavities contain the majority of fluid storage and provide fluid drainage to the fractures. Interporosity flows between the triple media are described using a pseudosteady-state concept and the system is characterized by interporosity transmissivity ratios and storativity ratio of each continuum. Pressure behavior is analyzed by examining the pressure drawdown curves, the derivative plots, and the effects of the characteristic parameters. Typical pressure responses from fractures, matrices, and cavities are represented by three semilog straight lines; the transitions by two troughs below the stabilization lines in the derivative plots. The analytical solution to the proposed model is further verified using a numerical simulation. The analytical model has also been applied to a published field-buildup well test and is able to match the pressure buildup data.

Parallel computing simulation of fluid flow in the unsaturated zone of Yucca Mountain, Nevada.

This paper presents the application of parallel computing techniques to large-scale modeling of fluid flow in the unsaturated zone (UZ) at Yucca Mountain, Nevada. In this study, parallel computing techniques, as implemented into the TOUGH2 code, are applied in large-scale numerical simulations on a distributed-memory parallel computer. The modeling study has been conducted using an over-1-million-cell three-dimensional numerical model, which incorporates a wide variety of field data for the highly heterogeneous fractured formation at Yucca Mountain. The objective of this study is to analyze the impact of various surface infiltration scenarios (under current and possible future climates) on flow through the UZ system, using various hydrogeological conceptual models with refined grids. The results indicate that the 1-million-cell models produce better resolution results and reveal some flow patterns that cannot be obtained using coarse-grid modeling models.

Modeling thermal-hydrological response of the unsaturated zone at Yucca Mountain, Nevada, to thermal load at a potential repository.

This paper presents a numerical study on the response of the unsaturated zone (UZ) system of Yucca Mountain to heat generated from decaying radioactive wastes emplaced at the proposed repository. The modeling study is based on the current thermal-hydrological (TH) mountain-scale model, which uses a locally refined 2D north-south cross-section and dual-permeability numerical approach. The model provides a prediction of the mountain-scale TH response under the thermal-load scenario of 1.45 kW/m, while accounting for future climatic changes and the effects of drift ventilation. The TH simulation results show that ventilation of the repository drifts has a large impact on thermal-hydrologic regimes and moisture-flow conditions at the repository. In both cases, with and without ventilation, the TH model predicts dry or reduced liquid saturation near the drifts for over 1000 years, during which liquid flux through the drifts is reduced to either zero or less than the ambient flux. Without ventilation, the model predicts higher temperatures at the repository, but no major moisture redistribution in the UZ except in the areas very near the heated drifts.

Characterization of flow and transport processes within the unsaturated zone of Yucca Mountain, Nevada, under current and future climates.

This paper presents a large-scale modeling study characterizing fluid flow and tracer transport in the unsaturated zone of Yucca Mountain, Nevada, a potential repository site for storing high-level radioactive waste. The study has been conducted using a three-dimensional numerical model, which incorporates a wide variety of field data and takes into account the coupled processes of flow and transport in the highly heterogeneous, unsaturated fractured porous rock. The modeling approach is based on a dual-continuum formulation of coupled multiphase fluid and tracer transport through fractured porous rock. Various scenarios of current and future climate conditions and their effects on the unsaturated zone are evaluated to aid in the assessment of the proposed repository's system performance using different conceptual models. These models are calibrated against field-measured data. Model-predicted flow and transport processes under current and future climates are discussed.